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Optical interferometers comprising birefringent-crystal beam displacers, wave plates, and phase shifters 
serve as stable devices for simulating quantum information processes such as heralded coined quantum 
walks. Quantum walks are important for quantum algorithms, universal quantum computing circuits, 
quantum transport in complex systems, and demonstrating intriguing nonlinear dynamical quantum 
phenomena. We introduce fully controllable polarization-independent phase shifters in optical pathes in 
order to realize site-dependent phase defects. The effectiveness of our interferometer is demonstrated 
through realizing single-photon quantum-walk dynamics in one dimension. By applying site-dependent 
phase defects, the translational symmetry of an ideal standard quantum walk is broken resulting in 
localization effect in a quantum walk architecture. The walk is realized for different site-dependent phase 
defects and coin settings, indicating the strength of localization signature depends on the level of phase due 
to site-dependent phase defects and coin settings and opening the way for the implementation of a 
quantum-walk-based algorithm. 

Quantum walks (QWs) 1 " 3 are the quantum mechanical analog of classical random walks (RWs), and hence 
can be used to develop quantum algorithms 4 " 7 , emerge as an alternative to the standard circuit model for 
quantum computing 8 " 10 , and represent one of the most promising resources for the simulation of 
physical system and important phenomena such as energy transport in photosynthesis 1112 , quantum chaos 13 " 16 , 
Anderson localization 17 " 23 , and topological phases 24 . 

A standard model of a one-dimension (ID) discrete-time QW consists of a quantum walker carrying a 
quantum coin. The walker goes back and forth along a line and the direction at each step depends on the result 
of a coin flip, which can be implemented by an arbitrary unitary operation in SU(2) following by a conditional 
position shift operation. The position variance of the walker a 2 = (x 2 ) — (x) 2 is linear on the number of the steps 
for RWs and quadratic for QWs. The position distribution of a standard QW P(x) shows a ballistic diffusion and 
that of RWs diffusive spreading. Furthermore if the static disorder is introduced in the dynamics of QWs, by 
changing the interference pattern localization effect can be observed in a QW architecture— spreading more 
slowly than RWs 19 . 

Experimental QWs began in 1999 and were performed with the frequency space of an optical resonator 2 . 
Several alternative realizations were quickly afterwards based on energy levels in nuclear magnetic resonance 25 , 
phase and position space of trapped ions 26 " 28 and trapped neutral atoms 29 , photons in beam splitter array, in fiber 
loop and in waveguide structures 3,16 ' 30 " 36 . 

In this work, we report on the implementation of a discrete QW with site- dependent single-point phase defects 
(SPPDs), and implement methods suggested by Wojcik et al 19 . We investigate the evolution of single-photons 
moving in a discrete environment presenting SPPDs, casting light on the natural feature of the physical descrip- 
tion of QWs such as general properties of diffusion modified by quantum or interference effects. Compared to the 
previous optical experiments on standard QWs (e.g. without static disorder) 3 ' 32 " 36 , we use phase shifters (PSs) for 
realization of the site-dependent phases acquired by the walker at the certain position. Compared to the experi- 
mental realization of QWs with time-dependent phase defects 23 , our experiment on QWs with site- dependent 
phases is more close to nature and can be used to investigate localization effect on low- dimensional structure, 
which would be interesting in research on properties of low-dimensional materials. 

Results 

The QW eliminates random evolution by trading the coin for a quantum two -level system, which, in our case, is 
the polarization state of a single photon: horizontal (H) and vertical (V). Furthermore the QW employs unitary 
dynamics by which the walker's position is entangled with the coin state. 



SCIENTIFIC REPORTS | 4:4825 | DOI: 1 0.1 038/srep04825 



1 



trigger 




1 



single photon 



BBO 




PBS 



HWP 



QWP 



PS 



OC 



BD 



APD 



Figure 1 | Experimental scheme for 10-step QW with site-dependent phase function. Single-photons created via type-I SPDC are injected to the optical 
network. Arbitrary initial coin states are prepared by a polarizing beam splitter (PBS), half- wave plates (HWPs) and quarter-wave plates (QWPs). 
Phase shifters(PSs) are placed in the corresponding spatial modes and the optical compensators (OCs) are used to compensate the temporal delay caused 
by PSs. Coincident detection of photons at avalanche photo-diodes (APDs, 7 ns time window) predicts a successful run of the QW. 



We use the array of beam displacers (BDs) as interferometer net- 
work similar to the setup in the previous work 16 . By taking advantage 
of the intrinsically stable interferometers, our approach is robust and 
able to control both coin and walker at each step. Benefiting from the 
fully controllable implementation, we experimentally study the 
impact of the SPPD and coin bias on the localization effect in a 
QW architecture and the experimental results agree with the theor- 
etical predictions. Compared to the previous experimental results 
which only simulated localization effect by trapping the walker in 
the original position x = 0, we experimentally localize the single- 
photons in different positions. 

In our experiment, we are able to achieve 10 QW- steps, surpassing 
the previous limit of 7 QW-steps in a similar interferometric setup 3 . 
The challenge of our experiment is to realize specific polarizing- 
independent phase on each site via microscope slides with precise 
effective thickness as PSs and to keep high interference visibility for 
each step even with phase defect. By introducing controllable PSs in 
paths of the interferometers, We have managed to create these ver- 
satile interferometer networks which can be used in many other 
fields. 

The coin state is encoded in the polarization \H) and |V) of the 
input single-photon which are generated by type-I spontaneous 
parametric downconversion (SPDC). In the basis {|H), | V)}, the coin 
operation for each step is given by 



C(B) = 



cos 20 sin 20 



sin 20 



- cos 20 



(i) 



with 0 G (0°, 45°), consisting of a polarization rotation, which is 
realized with a half-wave plate (HWP) setting. For example, the 
Hadamard operator is realized with a HWP set to 0 = 22.5°. 

The walker's positions are represented by longitudinal spatial 
modes which are implemented by birefringent calcite BDs. The 
optical axis of each BD is cut so that vertically polarized light is 
directly transmitted and horizontal light undergoes a lateral displace- 
ment into a neighboring mode which interferes with the vertical light 
in the same mode. Each pair of BDs forms an interferometer. 

The unitary operator for each step of a QW with a SPPD 



U= (S^®|H)(H|+SJ B ®|V)(V|) (fi®C(0)) 



(2) 



with the position shift operator % n \x) =e i ^ x ~ n) \x-\)[s^ n \x) = 

e x§b{x-n) |x+ 1)) on the modes manipulates the wavepacket to prop- 
agate according to the polarization of the photons. Here b(x) is the 



Dirac delta function which satisfies (5(0) = 1 and 5(x) = 0 for x 0. 
The translational symmetry of an ideal standard QW without SPPD 
is now broken by modifying the phase of the walker on each site, 
which can be realized by simply introducing PSs in the specific inter- 
ferometer arms. PSs are placed in the certain modes x = n. By 
adjusting the relative angle between the PS and the following BD 
the effective thickness of the PS changes and the specific phase (j) 
can be realized. 

The first 10 steps of the QW with SPPD (j) applied in the original 
position x = 0 are shown in Fig. 1 in detailed. The longitudinal spatial 
modes after the 1st step are recombined interferometrically at the 
2nd step. In our experiment, we attain interference visibility of 0.998 
(extinction ratio 1000 : 1) for each step, i.e., for each pair of sequential 
beam displacers. The photons emerge in the N + 1 spatial modes at 
the output of the Nth step and are subsequently detected by ava- 
lanche photodiodes (APDs). The probabilities are obtained by nor- 
malizing photon counts on each site to total number of photon 
counts for the respective step. The measured probability distribu- 
tions for 1 to 10 steps of a Hadamard QW with SPPD (j) = 180° and 
the antisymmetric initial coin state (\H) — i\V)) j \/2 are shown in 

Fig. 2a. An average distance d= \P exv (x) - P th (x)\ is 0.046 

ensuring a good agreement between the measured probabilities 
and theoretic predictions after 10 steps. The walker state after 4 steps 
clearly shows the characteristic shape of a localization distribution: a 
pronounced peak of the probability 0.615 ± 0.011 in the original 
position x = 0 and the low probabilities in the side positions. In 
contrast to the ideal standard Hadamard QW the expansion of the 
wavepacket is highly suppressed and the probability of the walker 
returning to the original position is enhanced strictly and displays the 
signature of the localization effect. 

In Figs. 2b and 2c, the position variance and recurrence probabil- 
ity, i.e., the probability of the walker returning to the original posi- 
tion, show the spread of the localized QW is much slower than that of 
the standard QW and the walker is trapped in the original position 
with high probability after the 4th step. While in the case of the 10- 
step standard Hadamard QW without SPPD the variance is given by 
Oq = 29.951, a lower variance occurs in the RW case with o\ = 10. 
Our measured value a 2 = 9.586 ± 0.463 agrees well with the theor- 
etical prediction 9.547 and shows an even slower spread than RW. 
The presented error bars include only statistical errors, calculated 
from the standard deviations of the values calculated by the Monte 
Carlo method. The measured recurrence probability exhibits the 
localization effect of QW with SPPD after the 4th step. Compared 
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Figure 2 | Localization effect in a QW architecture, (a) Experimental position distributions for successive steps of the Hadamard (C(22.5 0 )) QW with the 
SPPD (j) = 180° in the original position x = 0 and antisymmetric initial coin state up to 10 steps, (b) Measured position variance and (c) recurrence 
probability of localized QW for 1 to 10 steps, compared to theoretical predictions (solid lines). Error bars indicate the statistical uncertainty. Inset shows 
comparison of theoretical prediction of the recurrence probability of localized Hadamard QW with the SPPD 0 = 180° and standard QW with 0 = 0. 



to the walker in a standard QW with SPPD showing a ballistic beha- 
vior, it is always trapped in the position x = 0 with high probability 
about 0.64 after 4 steps shown in Fig. 2c. 

Our second experimental result highlights the full control of the 
implementation of the QW. In Fig. 3, we show the impact of single- 
point phase </> El [0°, 180°] on the localization effect. In this case, we 
change the effective thickness of the PSs to realize different SPPDs in 
the original position. Fig. 3a shows the position distribution of the 
10-step Hadamard QW changes as a function of the single-point 
phase (/). At the 10th step the recurrence probability Pio(0) and the 
position variance a 1 as functions of <j> are shown in Fig. 3b. For the 
antisymmetric initial coin state, the recurrence probability Pi O (0) 
increases with (j> in the range [0°, 135°] and decreases in the range 
(135°, 180°] . The maximal probability is achieved (measured as 0.660 
± 0.012 and numerical simulated as 0.667) with </> = 135°. The 
localization effect occurs in the range [45°, 180°], which agrees with 
the analytic result. Thus with fixed coin toss and initial state, whether 
or not the localization effect can be observed depends on the choices 
of single-point phase applied in the original position. 

The dependence of the localization effect on the phase </> can be 
explained 19 by the overlap between the localized eigenstates of the 
single- step unitary operator U and the initial state of the walker + 
coin system shown in Fig. 4a. The eigenstates and the corresponding 
eigenvalues can be found by solving a set of recurrence equations 
generated by the unitary operator U. The number of the localized 
eigenstates of U depends on (j). In the range <j> £ (0°, 45°) and 
<j> e (135°, 180°) there are two localized eigenstates. Whereas, in 
the range <j> G (45°, 135°) there are four such states. With the initial 



state \0)®(\H)-i\V)) j y/2, the overlap increases from 0 to 0.828 

with § G (0°, 135°) and decrease to 0.8 with <j> e (135°, 180°). The 
localization effect occurs in the range <j> £ [45°, 180°] and becomes 
most notable with </> = 135°. 

The next experimental result shown in Fig. 5 highlights the flex- 
ibility of our implementation with respect to the easy adjustability of 
the coin bias. Now we study the impact of the different coin biases on 
the localized QWs. Figure 5a shows the measured position distribu- 
tions for the 10-step localized QW with antisymmetric initial coin 
state, SPPD </> = 180° in the original position x = 0 and different coin 
biases 0 = 9°, 18°, 22.5°, 30° realized via different HWP settings. The 
walker is trapped in the original position with the proper choice of 0, 
which can also be observed from the position variances and recur- 
rence probabilities in Fig. 5b. With the angle of HWP 9 increasing, 
the mode of behaviour of the walker + coin system is transmitted 
from the diabatic (diabatic transition probability D = cos 2 20 ~ 1) to 
adiabatic limit (D^Cl). Harmin predicted that suppression of dif- 
fusion occurs in the form of almost perfect recurrences of the initial 
level population in the adiabatic limit, while in the diabatic limit, the 
recurrence will vanish 37 . Our experiment result agrees with the the- 
oretic predictions and shows that with the HWP angle of 0 = 9° and 
D = 0.905, in the diabatic limit the walker spreads widely and no 
recurrence occurs. Whereas, with the HWP angle increasing to 
0 = 30° (D = 0.095), the mode of the behaviour of the system is 
transmitted to the adiabatic limit, the diffusion is more suppressed 
and the recurrence probability increases. Thus the localization effect 
becomes more obvious in the adiabatic limit. This result can also be 




Figure 3 | The strength of localization signature due to the level of phase defects, (a) Experimental data of probability distributions of the 10-step 
Hadamard QW with antisymmetric initial coin state and various single-point phase (j) in the original position x = 0. The inset shows the probability 
distribution of the lOth-step Hadamard QW with SPPD 0 = 135°. Red and black bars show experimental data and theoretical predictions respectively, (b) 
Measured position variance of the 10-step Hadamard QW v.s. 0, compared to theoretical predictions (solid lines). The inset shows the measured 
recurrence probabilities Pio(0) after the 10 steps as a function of 0. 
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Figure 4 | Explanation of localization effect in a QW architecture, (a) Overlap between the localized eigenstates of the unitary operation U of the 
Hadamard QW with SPPD in x = 0 and the initial state of the walker + coin system as a function of the phase (j). (b) Overlap between the localized 
eigenstates of U of the QW with SPPD (j) = 180° in x = 0 and the initial state as a function of the coin bias 6. 



explained by the dependence of the overlap between the localized 
eigenstates and the initial state of the walker + coin system on the 
coin bias shown in Fig. 4b. With 0 increasing from 0° to 30°, the 
overlap increases monotonically from 0 to 0.974, which suggests the 
localization effect becomes more notable. 

Compared to the previous experimental results which only simu- 
lated Anderson localization of the QW by trapping the walker in the 
original position, we experimentally localize the single-photons in 
different positions, for example in x = 1. By inserting the PSs with 
proper effective thickness in the spatial mode x = 1 at each odd step, 
we can realize a localized QW with SPPD (j) = 180° in x = 1. In Fig. 6, 
the walker appears in x = 1 with highest probability after 9 steps. 
With different coin settings, the localization effect becomes more 
notable when the system is transmitted from diabatic to adiabatic 
limit. The tendency of the position variance and the probability of the 
walker back to x = 1 depends on the coin bias 0 in the same manner 
compared to the case of the localization in the original position. 

The performance of our setup is limited only by imperfections of 
the optical components such as nonplanar optical surfaces and the 
coherence length of single photons, resulting in errors and decoher- 
ence. The most significant source of systematic errors in our setup is 
the imperfect coherence visibility of the BD interferometer. A lim- 
itation for the maximal step number is given by the size of the clear 
aperture of BDs. For example, for 10-step QW, the effective diameter 
of the clear aperture of the BD with beam separation 3 mm needs to 
be larger than 30 mm. However this problem is not intrinsic to this 
implementation, since the BDs with large enough clear aperture and 
strictly planar surface can realize the large-step QW. 



Conclusion 

In summary, we implement a stable and efficient way to realize QWs 
embedded in a broader framework and show the phase defects can 
influence the evolution of wavepackets. The QW with SPPD has the 
single -photons localized in the certain position. Our experiment 
benefits from the high stability and full control of both coin and 
walker at each step and in each given position. The versatility of 
our setup allows for extensions, such as the realization of multi- 
particle QWs, in which richer choices of coin flipping and defects 
would help us to study the topology of arbitrary graphs and develop 
the applications such as quantum state transfer and energy trans- 
portation problems. The localization can also be used to filter and to 
trap particles, which would find applications in quantum algorithms 
and quantum state engineering. Our results show a new realm of QW 
phenomena and our new interferometer phase-shift control provides 
a valuable new tool for exploring QWs with various potentials. 
Furthermore, the versatile interferometer network with high inter- 
ference visibility can be used in many other fields of optical quantum 
information processing. 

Methods 

Photon-pair generation. The non-degenerated polarization degenerate photon pairs 
generated via type-I SPDC in two 0.5 mm-thick nonlinear-jS-barium-borate 
(BBO) crystals cut at 29.41°, pumped by a 400.8 nm CW diode laser (LBX-405-100, 
Oxxius) with up to 100 mW of power. For ID QWs, by triggering on one photon, the 
other at wavelength 800 nm is prepared into a single-photon state. A polarizing 
beam splitter (PBS) following by waveplates allow generation of any polarized state of 
single photon (e.g. any initial coin state). Interference filters determine the photon 
bandwidth 3 nm and then individual downconverted photons are steered into the 
optical modes of the linear-optical network formed by a series of birefringent calcite 




Figure 5 | Influence on the localization effect in a QW architecture due to different coin settings, (a) Experimental data of probability distributions of 
the 1 0-step QW with SPPD (j) = 1 80° in the original position x = 0, antisymmetric initial coin state and various coin bias 6. The inset shows the probability 
distribution of the lOth-step QW with the coin bias 6 = 30°. Red and black bars show experimental data and theoretical predictions respectively, (b) 
Measured position variance of the localized QW with antisymmetric initial coin state for 1 to 10 steps, with respective theoretical simulation (solid lines). 
The inset shows the measured probabilities of the walker returning to the original position Pi O (0) at the 10th step with the antisymmetric initial coin state 
as a function of the coin bias 6. Some of the statistical error bars are smaller than the symbol size. 
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Figure 6 | A 9-step QW with SPPD at x = 1. (a) Experimental data of probability distributions of the 9-step QW with antisymmetric initial coin state, 
SPPD (/) = 180° in the position x = 1 and various coin bias 6. The inset shows the probability distribution of the 9th-step QW with the coin flipping C(30°). 
Red and black bars show experimental data and theoretical predictions respectively, (b) Measured position variance of the 9-step QW v.s. the coin bias 6 y 
with respect to the theoretical predictions (solid line). The inset shows the measured probabilities of the walker returning to the x = 1 position P 9 (l). 



BDs, HWPs and PSs. Output photons are detected using APDs (SPCM-AQRH-14- 
FC) with dark count rate of <100 s" 1 whose coincident signals— monitored using a 
commercially available counting logic (id800-TDC)— are used to post-select two 
single-photon events. The total coincident counts are about 300 s" 1 (the coincident 
counts are collected over 60 s). The probability of creating more than one photon pair 
is less than 10" 4 and can be neglected. 

Beam displacer for interferometry. The spatial mode is implemented by a 
birefringent calcite BD with length 28.165 mm and clear aperture 33 mm X 15 mm. 
The optical axis of each BD is cut so that vertically polarized light is directly 
transmitted and horizontal light undergoes a 3 mm lateral displacement into a 
neighboring mode which interferes with the vertical light in the same mode. Each pair 
of BDs forms an interferometer. Only odd (even) sites of the walker are labeled at each 
odd (even) step, since the probabilities of the walker appearing on the other sites are 
zero. 

Ten BDs are placed in sequence and need to have their optical axes mutually 
aligned. Coalignment ensures that beams split by one BD in the sequence yield 
maximum interference visibility after passing through a HWP and the next BD in the 
sequence. 

Output photons from the interferometric network are coupled into a single-mode 
optical fibre and subsequently detected by a single-photon APD in coincidence with 
the trigger photon. We characterize the quality of the experimental QW by its 1-norm 

distance 3 from the simulated QW according to d= ^ |P exp (%) -P th (x) |. The 

distance increases monotonically with each step number due some lack of of relative 
phase control between the multiple interferometers. This limited control is due to 
nonplanar optical surfaces. 

Phase shifters. We introduce microscope slides into the interferometric network after 
first completing the alignment described above and ensuring maximum interference 
and small distance between simulated and empirical walker distributions. The 
microscope slides are inserted and aligned to recover the case of zero phase shift 0 = 0 
for all locations. This microscope- slide alignment corresponds to each slide being in 
the plane perpendicular to the beams in each path. 

After achieving this alignment, each microscopic slide can be adjusted to an 
effective thickness in order impart a controllable PS independent from all other PSs. 
This effective thickness is achieved by rotating the slide out of the plane perpendicular 
to the beams. Tilting the slide is a viable alternative but is not as stable for the times 
required to gather the data. 
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